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ABSTRACT 


Progress in computers during the past two decades has sparked the development 
of many useful high frequency (HF) ionospheric propagation prediction codes. The high- 
latitude (polar) ionosphere still remains as the most difficult propagation region to 
predict. 

A "Noncentric" database of collected high-latitude signal and noise measurements 
was obtained during 1988 and 1989 by the University of Leicester, U.K.. 

The Advanced Prophet 4.0 HF Ionospheric Propagation Prediction Code was 
exercised and compared to a portion of the "Noncentric" database, for a transmitter at 
Clyde River, Canada and a receiver at Leicester, U.K. The Prophet predictions were 
better during winter months than during summer months. Overall, 70% of the signal 


strength data from Prophet 4.0 were between -20 dB and +20 dB error, compared to the 


measured data. 
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I. INTRODUCTION 


A. OVERVIEW 

High frequency (HF) (2-30 MHz) electromagnetic waves have 
been, and still are, extensively used to meet the plethora of 
military and civilian communications requirements. 

Although in the 70's HF fell into disfavor because 
communication satellites emerged as an alternative to HF, in 
the 80's HF regained its vital role as aé flexible, 
inexpensive, and resilient ccmmunications capability. HF 
offers low cost equipment, low power requirements, and 
adequate bandwidths. 

From the military point of view, reliable communications 
between the decision making centers and units in the air, sea, 
and land are essential. 

For large military organizations like N.A.T.O. long- 
distance communications are extremely important since the area 
of operations extends from the high Arctic to the equator and 
from the Mediterranean Sea to the Pacific Ocean. 

Long-distance communications at HF is achieved via 
skywave propagation. Skywaves are launched at a range of 
angles and travel upwards to be reflected from the ionosphere. 


The ionosphere varies diurnally, seasonally, and with solar 


activity in observable patterns, but is still quite 
unpredictable and variable for short-term applications. 

Of particular interest is the high-latitude ionosphere 
which is controlled to a great extent by the earth's magnetic 
field. The variability of high-latitude signals during 
disturbances can be attributed to the variability in shape and 
density of high-latitude ionosphere. 

For many years high~latitude HF propagation had been an 
enigma primarily due to insufficient analyzable data. During 
the last decade, however, many experiments have helped to 
explain and permitted analysis of some high-latitude 
lonospheric phenomena. 

In Appendix A, the kasic concepts of skywave propagation 
and the effects of the high-latitude ionosphere on HF 


communications are briefly presented. 


B. PROPAGATION MODELS 

The remarkable progress in computer systems over the past 
two decades has made it practical to develop models of the 
ionosphere with the objective of making propagation 
predictions. 

There are two kinds of ionospheric models, empirical and 
mathematical. Empirical models are based on data compilation; 
conversely, mathematical models are based on thevretical 


physical descriptions of the ionosphere. 


Propagation models are constructed by adding ray-tracing 
routines, attenuation models, and noise models to an 
ionospheric model in order to predict the characteristics of 
propagation (field strength, maximum usable frequency, 
propagation mode etc.). The ioncspheric model depends on an 
empirical database and the propagation models, therefore, omit 
day-to-day variability and tend to be less accurate in 
equatorial and polar regions. 

A large number of prediction codes are available such as: 
PROPHET, IONCAP, WOMAP, AMBCOM, HFMUFES, APPLAB, CCIR252 etc. 
This thesis is concerned with the Advanced PROPHET 4.0 [1] 
propagation prediction code and the accuracy of its 
predictions for polar paths. In Chapter II, the capabilities 


of Advanced PROPHET are described. 


C. THE CLYDE RIVER-LEICESTER SHORTWAVE LINK 

Recently, scientists from the University of Leicester, 
U.K., conducted a series of signal and noise measurements on 
a shortwave link from Clyde River (CANADA) to Leicester, U.K., 
(Figure 1). A “Noncentric" database thus created may be the 
largest collection of polar signal and noise data to date. The 
"Noncentric" term is used only for identification of data 
source, since these measurements were a subset of data from 
Project "Noncentric". 

The "Noncentric" database was not constructed for the 


evaluation of PC-based propagation prediction codes. It is not 





Figure 1. The geographical location of the Leicester-Clyde 
River shortwave link. 


applicable for incoming field strength comparisons, but for 
signal-to-noise ratio comparisons at receiver input terminals. 
The reasons for the database not providing incoming field 
strength are twofold: First, the antenna. radiation 
characteristics at the transmitter and at the receiver were 
not established and second, the RF signal distribution system 
characteristics between the receive antenna terminals and the 
input to the receiver were never determined, and if they had 
been, it is highly probable that they would have varied 
substantially with time. 


The "Noncentric" database covers 1988 and 1989 but only 
data of “site O" (summer 1988 and winter 1989) were available 
and utilized in this study. 

The receiver antenna, located at Leicester, is a sloping- 
vee antenna which consists of two terminated sloping wire 
antennas with a 90° interior angle. Chapter III describes how 
the antenna was simulated and analyzed using the Numerical 
Electromagnetics Code (NEC) [2]. The calculated radiation 
patterns were installed in Advanced PROPHET for the receiving 
site. 

The transmitter antennas, located at Clyde River, were a 
loaded vertical monopole for the winter 1989 campaign, and a 
"Butternut" trap-vertical monopole for the summer 1988 
campaign. Their radiation patterns were obtained from Ref. 3 
and were also installed in Advanced PROPHET for the 


transmitter site. 


D. ANALYSIS OF PREDICTIONS 

In Chapter IV, Advanced PROPHET's predictions of signal- 
to-noise ratio were compared with the University of 
Leicester's measured signal-to-noise data, by computing the 
prediction error, e, of each campaign, where e = (S/N) pao — 
(S/N) e¢- Finally, the predictions of the two campaigns were 


compared to show the seasonal variation. 


II. DESCRIPTION OF ADVANCED PROPHET 4.0 


Advanced PROPHET 4.0, developed at NOSC in 1989, is a 
successor to Advanced PROPHET 3.2. Advanced PROPHET is a 
collection of computer simulation models designed to support 
tactical use of the HF band. 

It can determine the existence of an HF skywave channel 
between a transmitter and a receiver located anywhere in the 
world (1]. The modeling approach in PROPHET uses empirical 
models which have been fitted to experimental datz. That is, 
the models were derived from an extensive data set of 
measurements taken from various sources rather than by 
applying a detailed physical model. 

Advanced PROPHET consists of several computer algorithms 
capable of solving a great variety of HF problems. The 
standard PROPHET skywave calculation method involves three 
important algorithms: "MUFLUF", "MUF85" and "DAMBLT". 

The main objective of any HF prediction program is to 
determine accurately the skywave propagation bandwidth between 
the LUF and the MUF for any set of conditions. This is 
achieved by the “MUFLUF" algorithm which computs the path 
geometry between the transmitter and receiver. The MUF is 


estimated at the path midpoint using the "MUF85" algorithm. 


"MUF85" (also called MINIMUF) is a curve fit of first order 





exponential functions to a large database of HF sounder 
measurements. It provides a good approximation of MUF between 
any two geographical points. The "classical MUF" is scaled by 
interpolating a data table to determine the operational MUF 
and also the frequency of optimum transmission (FOT). The FOT 
is defined as a frequency with 90% reliability. A comparable 
algorithm has been developed to estimate the LUF. The LUF is 
determined in part by the solar insolation along a great 
circle path between two points. The algorithm for the 
estimation of LUF has been improved in version 4.0. 

Once the operational MUF and LUF have been established, 
the signal strength is determined by the "DAMBLT" algorithn. 
Constructed by T. Damboldt of the German Bundesposte, “DAMBLT" 
is an empirical fit of signal strength. This empirical fit 
depends upon the ratio of the operating frequency to the MUF, 
the ratio of the LUF to the frequency, and the path length. 
PROPHET 4.0 is also equipped with an improved algorithm for 
field strength predictions in polar regions. 

The "MUFLUF", “MUF85", and "DAMBLT" models comprise the 
basis for PROPHET's predictions. A number of other empirical 
models offer more sophisticated capabilities: (1) the ray- 
tracing program that computes the propagation modes and the 
actual path traversed by an HF signal between the transmitter 
and receiver; (2) an empirical model that locates the auroral 


oval with mapping options; and (3) an empirical model for the 


scintillation of downlink signals from a geostationary 
satellite. 

A model for the expected noise for various environments 
includes contributions from atmospheric, man-made, and 
galactic noise sources. Influences by jammer sites are also 
included. 

A large variety of input parameters form the different 
scenarios of HF problems. Parameters of great importance 
include: (1) antenna radiation patterns, (2) modulation 
parameters such as modulation type, bandwidth, and required 
signal-to-noise ratio, (3) geophysical parameters such as 
solar sunspot activity, magnetic and ionospheric storms, and 
(4) environmental parameters such as the type of terrain and 
cover for groundwave calculations. (Groundwave calculations 
with PROPHET are not recommended though, until existing errors 
are corrected). 

Advanced PROPHET products are utilized to guide operators 
of HF communication equipment to improve the effectiveness of 
the HF channel, and to reduce the probability of intercept and 
the vulnerability to a jammer threat. 


III. THE SLOPING-VEE ANTENNA MODEL 


A. INTRODUCTION 

Long-wire antennas are usually several wavelengths long 
and are often used in the HF band (2 to 30 MHz). The best 
known types are: the horizontal-vee antenna, the horizontal- 
rhombic antenna, the sloping-rhombic antenna, the vertical-vee 
antenna, the vertical inverted-vee antenna, the sloping-vee 
antenna, and the single horizontal wire antenna. 

Long-wire antennas may be characterized as resonant 
antennas (standing-wave) and traveling-wave antennas. In the 
case of resonant antennas, the current on the wire is a 
standing wave with sinusoidal variation. Resonant antennas 
usually operate satisfactorily only at a particular frequency 
and harmonics of this frequency. The input impedance is very 
dependent on frequency, resulting in narrow-band operation. 

For the traveling-wave case, a long-wire antenna is 
operated as a nonresonant antenna by virtue of termination of 
the far end of the wires using a load resistor which matches 
the characteristic impedance of the antenna when it is viewed 
as a transmission line. Traveling-wave antennas have wide 
bandwidth and present a useable impedance over a wide range of 


frequencies. 


B. THE MODELING PROCESS 

The sloping-vee antenna, used at the University of 
Leicester U.K. as a receiving antenna in the Clyde River- 
Leicester link, consists of two sloping straight wire antennas 
arranged with an interior angle of 90°. Figure 2 illustrates 
the configuration and dimensions of the sloping-vee antenna. 
The two long arms of the antenna are terminated in 600 Ohm 
loads and a 1 meter ground stake, producing traveling-wave 
operation. 

The sloping-vee was modeled via NEC [3]. NEC was developed 
at the Lawrence Livermore Laboratory at Livermore California 
and is a user-oriented computer code for analysis of the 
electromagnetic response of antennas and other metal 
structures. It is a versatile code with wire antenna modeling 
capabilities that, include transmission lines, networks, 
loading and ground effects. 

A typical input data set is in Appendix C. Modeling a wire 
structure in NEC involves both geometrical and electrical 
factors. The main electrical consideration is the wire segment 
length relative to the wavelength. 

Forty-one segments were used for each of the 100 m long 
arms of the sloping-vee antenna. For the frequency band of 
interest, 3-20 MHz (15m<A<100m), the length of each segment 
varies from 0.024 to 0.16A over the band. These limits are 
within NEC modeling guidelines. The segment size determines 


the resolution for the current on the model since current is 
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ANTENNA LEGS (100m) 


0.6 m 
LOAD 600 2 
GROUND STAKES Im 
“GOOD GROUND” 





Figure 2. Configuration of the sloping-vee antenna 


computed at the center of each segment. 

Parameters for "good ground" (o=30, €=0.01) were used to 
simulate the terrain of the sloping-vee. Due to the presence 
of finite ground, the Sommerfield/Norton method was utilized 
for accurate results. 

Two feed points and the voltage source option were used to 
specify the excitation of the structure. The two feed points, 


seen in Figure 2, are located 1.2 m from the beginning of each 


long arn. 


i1 


\" 
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C. ANALYSIS OF THE RESULTS 

Figures 3, 4, and 5, contain the elevation and azimuth 
radiation patterns of the sloping-vee for 3, 8, and 20 MHz. 
The patterns are as expected from antenna theory. 

Finally, Figure 6 illustrates the variation of the input 
impedance over the band. Note that the variation is rather 
small, and the values of input impedance are within the 


expected range. 
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Figure 5. Elevation and azimuth radiation patterns of the 
sloping-vee antenna at 20 MHz. 
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VARIATION OF INPUT INPEDANCE 
WITH FREQUENCY 


Freq [MHz] 





Figure 6. Variation of the real and imaginary part of the 
input impedance of the sloping-vee antenna with the 
frequency. 








IV. COMPARISON OF PROPHET'S PREDICTIONS WITH MEASUREMENTS 


A. THE “NONCENTRIC" DATABASE 

The "Noncentric" database developed at the University of 
Leicester, U.K., consists of more than 80,000 signal and noise 
level data covering the two-year period of 1988-1989, for the 
Clyde River-Leicester path. 

The data was first characterized by "campaigns". The two- 
year period was divided into twelve campaigns of approximately 
25 days each from two consecutive months. The winter 1989 
campaign, for example, consists of data obtained during the 
time period from 1-17-89 to 2-12-89 (25 days); and, similarly, 
the summer 1988 campaign consists of data obtained from 7-17- 
88 to 8-12-88 (also 25 days). 

In order to speed-up the analysis of "Noncentric" data, 
two tests were developed [4] to detect the presence of the 
"Noncentric" signal. These are the Spread-Index test and the 
Call-Sign test. 

The incoming signals at the receiver were processed using 
a 1,000 point fast Fourier transform (FFT) to calculate the 
doppler spreading of the received signal. 

The Spread-Index test is based on ionospherically induced 
doppler spreading of a CW transmission. The spectrum of the 


received signal is normalized by setting the peak amplitude to 
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unity. The entire spectrum interval of the CW transmission is 
50 Hz (from -25 Hz to +25 Hz) as can be seen in Figure 7a. The 
mean noise level of the first quarter of the spectrum (-25 Hz 
to -12.5 Hz) is assumed to represent the mean noise level of 
the whole spectrum. Thus, it is subtracted to yield a spectrum 
as depicted in Figure 7b. The Spread-Index is the area under 
the remaining signal which is proportional to the doppler 
spreading of the signal. The Spread-Index, in conjunction with 
the signal level measurements, provides the first recognition 
test as follows: if the Spread-Index is large the received 
signal passes the test, but, if the Spread- Index is small, 
the signal passes the test only if its signal level is high. 

The Call-Sign test is based on the recognition of the call 
sign, CZB, transmitted three times in each two-minute 
sequence. A typical call-sign spectrum is illustrated in 
Figure 8. The signal recognition routine searches within the 
frequency interval from -5 Hz to +20 Hz in order to find a 
pattern similar to that of Figure 8 (a main signal peak with 
two side peaks). If this pattern is found, the signal passes 
the Call-Sign test. 

The data, after undergoing the two recognition tests, were 
categorized into the following groups: 

(1) Signals passing both tests are identified as group 

SIG. 


(2) Signals failing both tests are identified as group 


NOS. 





eeoerervraneasae sense esevna eaeeoeceeove 





Pigure 7. Frequency spectrum of a CW signal; (a) with 
noise, (b) with the noise removed [From Ref. 4]. 


(3) Signals passing the Call-Sign test only are 
identified as group FCS. 

(4) Signals passing the Spread-Index test only are 
identified as group FSI. 

(5) Signals characterized as interference are 


identified as group INT. 
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The "Noncentric" data base is comprised of the sum of all 
the groups except that from the INT group (interference). 
Overall, however, the confidence in signals that passed both 
tests (SIG group) is much higher than that of the other 


groups. For this reason, in the present thesis, only the SIG 
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Figure 8. The frequency spectrum of the call sign portion 
of the transmission (Ref. 4]. 


groups of the winter 1989 and summer 1988 campaigns were 


considered. 





B. PROPHET'S SCENARIO PARAMETERS 

In order to satisfactorily utilize the Advanced PROPHET's 
capabilities, an accurate creation of the appropriate scenario 
was necessary. 

The antenna parameters were of great importance. Radiation 
patterns of three different antennas were used in Advanced 
PROPHET; a loaded vertical monopole, a "Butternut" trap- 
vertical monopole as transmitting antennas, and a sloping-vee 
antenna as a receive antenna. 

In Appendix B, Advanced PROPHET's requirements for the 
three antennas are described. These patterns are listed in 
terms of the total gain (total = vertical + horizontal 
polarization) over the entire HF band (2-30 MHz). 

Geomagnetic and solar scenario parameters were taken from 
Ref. (6],{7]}] and [(8],(9] for the winter 1989 and summer 1988 
campaigns, respectively. The daily average of the eight, 
three-hourly indices, K,, was used as the input value for the 
planetary index K,- Average values of sunspot number measured 
by two different agencies were used. The two agencies were the 
American Association of Variable Star Observers and the WDC-A 
for Solar-Terrestrial Physics. 

The modulation type, CW with 50 Hz bandwidth, was also 
part of the Advanced Prophet scenario. 


Finally, atmospheric noise and man-made noise vere 


selected to simulate noise environments. 








C. ANALYSIS OF PREDICTIONS 


1. Assumptions 

In order to develop a quantitative measure of the 
prediction accuracy of Advanced Prophet, the model predictions 
of the signal-to-noise ratio, (S/N),,,, were compared with the 
empirical signal-to-noise ratio, (S/N) eres from the 
*"Noncentric" database. 

The statistical analysis of the prediction error, 
e = (S/N) pepo - (S/N)ieiee WAS Computed using the frequency 
distribution of the predicted errors and the cumulative 
distribution of the absolute value of the prediction errors. 

Furthermore, a relative measure of the prediction 
accuracy was achieved by normalizing each set of prediction 
errors to zero statistical mean and recomputing the frequency 
distribution and the cumulative distribution of the prediction 
errors. In this case, the prediction error, e', is given by 
e'= (S/N) peg ~ (S/N)ie¢ — E(e) , where E(e) is the mean of e. 

The computation of signal-to-noise ratio, rather than 
the field strength, was necessary due to the limited 
characteristics of the "Noncentric" signal strength data, and 
the lack of information about the complexity of the receiver 
site. 

The Clyde River-Leicester path is approximately 3,810 
Km which implies the existence of multihop propagation modes. 


Figure 9 is an examble of the ray-tracing routine of Advanced 
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Figure 9. Propagating modes in the Leicester-Clyde River 
shortwave link, at 10 MHz, on 8-9-1988 at 21.01 U.T.. 


PROPHET which was utilized to yield the maximum number of 
propagating modes for 10 MHz on 8-9-1988 at 21.01 Universal 
Time. 

The winter 1989 and summer 1988 campaigns involve 16 
HF frequencies ranging from 3.185 to 23.169 MHz. Each 
frequency was transmitted at a pre-specified minute, hourly. 
For example, the 13.886 MHz frequency was transmitted at 22 


minutes past the hour, the 10.195 MHz at 43 minutes past the 








hour, the 20.3 MHz at 49 minutes past the hour and so on for 
every hour of the day and for every day of the campaign. 

Since the prediction values of Advanced Prophet refer 
only to integer frequencies and at integer hours of the day, 
manual interpolation was conducted. 

A total of 2,348 data points was used in the 
evaluation of Advanced Prophet 4.0 (approximately 60 % of the 
data from the SIG files of winter 1989 and summer 1988 
campaigns). 

2. The Winter 1989 Campaign 

The SIG file of the winter 1989 campaign consists of 
1,979 data and 13 frequencies. In the present analysis, 1,127 
data were utilized for 6.8 MHz, 6.9 MHz, 9.9 MHz, 13.9 MHz, 
17.5 MHz and 20.3 MHz. Advanced PROPHET predicted useable 
paths in 95.8% of the cases (1,080 out of 1,121 cases). 

In Figure 10, the frequency distribution of the error, 
e, is depicted. The median error is -7.7 aB, and the standard 
deviation of the error is 14.7 dB. In Figure 11, the frequency 
distribution of the prediction error, e', may be seen. The 
median error is now only 1.98 dB since the statistical mean 
is, by definition, zero. 

Figure 12 illustrates the cumulative distribution of 
the absolute value of the error, e, (solid curve), and the 
cumulative distribution of the error, e', (dashed curve). It 


is observed (for the case of e) that 82% of Prophet's data are 
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Figure 10. Frequency distribution of the prediction error, 
e, for the winter 1989 campaign. 
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Figure 11. Frequency distribution of prediction error, e', 
(mean=o), for the winter 1989 campaign. 
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Figure 12. Cumulative distribution of the absolute value 
of Prophet prediction errors for the winter 1989 campaign. 








between ~20 dB and +20 dB error while (for the case of e') 


90.5% of Prophet's data are between -20 dB and +20 QB error. 


3. The Summer 1988 Campaign 

The SIG file of the summer 1988 campaign consists of 
1948 data points and 14 frequencies. Now, 1,221 data were 
utilized, involving the same frequencies as in the winter 1989 
campaign plus the 10.1 MHz frequency. Advanced Prophet 
predicted useable paths in 1,080 out of 1,221 cases (88%). 

Figures 13 and 14 show the frequency distribution of 
the predicted errors, e, and e'. The standard deviation of the 
predicted error is 14.28 a@B, and the median error is -21 dB 
and -0.24 dB for the prediction errors, e, and e', 
respectively. 

Figure 15 illustrates the cumulative distribution of 
the absolute value of the prediction errors, e, and e'. For 
the case of e (solid curve), 56% of Prophet's data are between 
-20 dB and +20 dB error. On the other hand, for the case of et 
(dashed curve), 88% of Prophet's data are between -20 GB and 
+20 GB error. 

4. Seasonal Variation of Predictions 

In Figures 16 and 17, the frequency distribution of 
the prediction errors, e, and e' is shown respectively for the 
two campaigns. The peak value of the error frequency for the 


winter 1989 campaign is higher than the peak value of the 
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Figure 13. Frequency distribution of the Prophet predicted 
error, e, for the summer 1988 campaign. 
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Pigure 14. Frequency distribution of the Prophet 
prediction error, e', for the summer 1988 campaign. 
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Figure 15. Cumulative distribution of the absolute value 
of Prophet prediction errors, e, and e', (summer 1988). 
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summer 1988 campaign by 6% (26% and 20% respectively). In 


Figure 16, for prediction error, e, the error value 
corresponding to the peak value of error frequency, winter 
1989 campaign, is -1 dB, and, summer 1988 campaign, -24 dB. 
Similarly, in Figure 17 for the prediction error, e', the 
respective numbers are 8 dB for the winter 1989 campaign and - 
3 GB for the summer 1988 campaign. 

Figure 18 shows the cumulative distribution of the 
absolute value of the prediction error, e, for the two 
campaigns. It is seen that 82% of the Prophet data of the 
winter 1989 campaign and only 56% of the Prophet data of the 
summer 1988 campaign are between -20 dB and +20 dB error. 

Finally, in Figure 19 the cumulative distribution of 
the absolute value of the prediction error, e', is depicted 
for the two campaigns. Note that 90.5% of the winter 
predictions and 88% of the summer predictions are between -20 
@B and +20 dB error. 

From the above analysis, it is concluded that the 
accuracy of Prophet 4.0 predictions is better during winter 
months than during summer months for the “Noncentric™ data 
analyzed. 

5. Macroscopic View of All Data Analyzed 

Figure 20 show the frequency distribution of the 

predicted error, e. The standard deviation of the error is 


15.23 a@B, and the median error -14.5 dB. In Figure 21, the 
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Figure 16. Frequency distribution of the Prophet 
prediction error, e, for the two campaigns. 
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Figure 17. Frequency distribution of the Prophet 
prediction error, e', for the two campaigns. 
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Figure 18. Cumulative distribution of the absolute value 
of the prediction error, e, for the two campaigns. 
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Figure 19. Cumulative distribution of the absolute value 
of the Prophet prediction error, e', for the two 
campaigns. 
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frequency distribution of the predicted error, e', for all the 


data analyzed, may be seen. The median error is now 0.77 GB, 
and the standard deviation of the predicted error is 14.17 dB. 

Finally, Figure 22 illustrates the cumulative 
distribution of the absolute value of the predicted error, e, 
(solid curve), and the cumulative distribution of the 
predicted error, e', (dashed curve). It is seen (for the case 
of e) that 70% of the predicted data are between -20 d@B and 
+20 dB error while (for the case of e') 90% of Prophet's 


predictions are between -20 dB and +20 GB error. 
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Figure 20. Frequency distribution of the Prophet 
prediction error, e, for all the data analyzed. 
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Figure 21. Frequency distribution of the Prophet 


prediction error, e', for all the data analysed. 
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Figure 22. Cumulative distribution of the absolute value 
of Prophet prediction errors, e, and e'. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

The sloping-vee receiving antenna was simulated and 
analyzed using the Numerical Electromagnetics Code (NEC). The 
radiation patterns and input impedance were as expected from 
the theory of traveling-wave antennas. 

Advanced Prophet 4.0 predictions are, in general, more 
accurate in winter than in summer for polar paths. 

Considering the variability of the polar ionosphere, the 
accuracy of Advanced Prophet's predictions for the simulated 
Clyde River-Leicester shortwave link is reasonable. 

Further evaluation of propagation prediction programs for 
polar paths should be conducted. Some of the programs are: 

The propagation model AMBCOM is available at NPS. It has 
an algorithm for computing the electron density profile in the 
high-latitude ionosphere and a model for auroral absorption, 
and is currently under evaluation. 

A new version of IONCAP-PC for polar regions, ICEPAC, is 
also available and will be evaluated during 1991. 

The "Noncentric" database provides a long-needed source of 
polar data for the evaluation of the above prediction models. 
However, the "Noncentric" data are not measurements of the 


incoming field strength but of signal-to-noise at the receiver 


input terminals. This is a major "shortcoming" of the database 





which will affect the quality of the evaluation of prediction 
models. 

It is recommended that any future data collected must 
accurately relate received field strength to receiver input 
voltage. 

A database involving several types of polar paths would be 
better for the evaluation process. The areas of interest for 
the location of transmitter and the receiver are: the polar 
cap, the auroral oval, the auroral trough, or outside of the 
above three. Thus, a total of 16 different paths are possible 
and should be consider carefully in planning future evaluation 


of high-latitude prediction codes. 


42 


APPENDIX A ELEMENTS OF SHORTWAVE THEORY 


A. GROUND WAVE PROPAGATION 

The lower part of the HF band (2-5 MHz) can be used for 
communications by means of ground wave. Ground wave theory 
ignores the effects of the ionosphere and assumes that the 
propagation path is around a spherical earth with homogeneous 
parameters. Ground wave is affected primarily by the 
electrical characteristics of the ground (conductivity and 
dielectric constant), climate, vegetation and the polarization 
of the wave (vertically polarized fields are most suitable for 
ground wave propagation). A characteristic feature of ground 
wave is that it is independent of diurnal and seasonal 
variations. The ranges obtained by ground wave propagation are 
determined by the attenuation of fields and are strongly 
dependent on the operating frequency. (Lower frequencies 


propagate over greater distances). 


B. SKYWAVE PROPAGATION 

Skywave propagation is more important than ground wave 
since it supports very long distance paths which often reach 
several thousand kilometers. The transmitted wave utilizes 
certain properties of the ionosphere as it is reflected back 
to earth. 
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1. Structure of the Ionosphere 


The ionosphere, located up in the stratosphere, 
blankets the earth and consists of several electrically 
conducting layers that vary widely over the earth's surface. 
The formation and intensity of these ionized layers are 
primarily due to electromagnetic radiation from the sun 
extending over the ultra-violet and X-ray portions of the 
spectrum. Here sufficient quantities of electrons and positive 
ions reflect the radio waves. The ability of the various 
layers to reflect HF waves becomes stronger as the radiation 
from the sun increases; at the same time, though, the 
attenuation of the transmitted wave also increases. 

The solar radiation varies considerably with the 
geographic latitude, the time of day, the season and the 
sunspot activity. The immediate effect of this variability is 
a change in the possibility of propagation being supported for 
certain frequenciege. The ionosphere is divided into three 
regions, or layers, designated as D, E and F layers. 

a. The D-layer 

The D-region is located at an altitude of 60 to 90 
km. It has a comparatively low free- electron density of about 
10° e/m’, compared with a much higher molecular density of 
about 10” mol/m’. The D-layer is a day-time phenomenon. It 
appears after sunrise and attains its strongest ionization 


shortly after local noon. During afternoon and evening hours 


44 





the layer gradually breaks down and is neutralized after 


sunset. There is also a seasonal variation in the D-layer 
electron densities with maximum value in summer. 

D-layer ionization is not sufficient to reflect 
shortwaves therefore penetration occurs. The relatively high 
electron collision frequency in this region is the main cause 
for the absorption on the D-layer. This attenuation of 
shortwaves in the D-region during daytime hours determines 
achievable ranges and transmitter powers necessary for a 
transmission. 

-laye 

The E-layer spans the altitude range from 90 to 130 
Km above the earth. Maximum ionization takes place at 110 Kn. 
The E-layer has a pronounced diurnal variation with maximum 
electron density around noon. At night, only a small residual 
level of ionization remains and it has virtually ceased to 
have any effect in shortwave propagation. The distances 
attained via the E-layer for daytime propagation are less than 
2000 Kn. 

In addition to the regular E-layer, there is also 
the sporadic-E layer, E,, named after its transient and 
irregular nature. It is situated 120 Km above the earth and 
exhibits very high ion concentration. Its occurrence depends 
largely on the geographical latitude. Thus, it is essentially 


a daytime phenomenon in low latitudes with no significant 





seasonal variation. Conversely, in high latitudes it occurs 
only at nighttime. In middle latitudes the formation of E, is 
more frequent by day than by night and more in summer than in 
winter. Because of these variations, the study of sporadic-E 
is done by statistical means. The high ionization of E, 
enables frequencies well above 30 MHz to be transmitted (VHF 
band). 
cc. The F-layer 

The F-layer is very important for shortwave 
communications because of the large attainable distances. It 
is divided into the Fl and F2 sublayers. 

The Fi-layer occupies the altitude region from 150 
to 220 Km above the earth and exists only during the day. It 
is more pronounced in summer than in winter, and at high 
sunspot numbers. The Fl-layer occasionally reflects 
shortwaves, but, usually, they penetrate the Fl region and are 
reflected by the F2-layer. 

The F2-layer is located 225 to 450 Km above the 
earth and is the principal reflecting region for long distance 
shortwave communications. In contrast to the E and Fl-layers, 
the F2-layer does not disappear after sunset and a residual 
ionization remains. Thus, during the night, HF propagation is 
only possible because of the F2-layer. The height and 
ionization of the F2-layer vary diurnally and seasonally with 


maximum height in summer and minimum in winter. Figure 2 shows 
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Figure 23. Diurnal variation of electron density and the 
formation of ionospheric layers. From [Ref. 12: p. 21). 


a typical distribution of the electron density at noon and 


midnight hours for an undisturbed ionosphere. 


2. The Propagation of HF-Waves in the Ionosphere 
ord Ww s 
The earth's magnetic field deflects the free 
electrons in the upper atmosphere from their linear path in 
such a way as to cause them to move in the form of a spiral 
around the magnetic lines of sGrce:. 
The strength of the magnetic field B (Wb/m*), as 


well as the mass m (gr), and the charge of the electron e 
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(Cb), are the quantities that determine the angular velocity 


n,, of the electron's spiral movement, 
Ny = (e/m)B (1) 


From (1) the frequency is easily obtained (n, = 
6.282f, ) and is called gyrofrequency f,. In the ionosphere 
f, = 2.84x10"B and varies from 1.7 MHz at the magnetic poles 
to 0.7 MHz at the magnetic equator [10]. 

When an electromagnetic wave of frequency f enters 
the ionosphere, it is split into two components: the ordinary 
f,, and the extraordinary f, wave. These are given respectively 


by 


rh 
° 
ut 


[f(f£+£,))'? , (2) 


and 


rh 
by 
il} 


[f(f£ - f))'% . (3) 


The polarization of the ordinary and extraordinary 
waves is circular, moving in opposite directions. The 
variability of the ionosphere affects the two waves. Their 
reflections take place at different altitudes; hence, the 
refractive indices, the delay times and the attenuation are 
different for the two waves. The extraordinary wave is 


attenuated more than the ordinary wave. 
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The two reflected waves, after leaving the 
ionosphere, recombine to form the originally transmitted wave 
of frequency f. Due to the differing attenuation and delay 
times of the two waves, the polarization of the emitted wave 
now becomes elliptical. 

A c Frequenc FO and _ LUF 

Snell's law of refraction shows, that, if a skywave 
is incident upon the ionosphere at angle i, with respect to 
the normal, the height at which the reflection occurs has a 
refractive index, n,, equal to sin(i,). For vertical incidence 
i, = 0, which means that sin(i,) = 0, and thus the Petiection 
will occur at the height where n, = 0. It is known [10), that 


n, is given by 
n, = (1 - (81N, /f? ))'? , (4) 


where N, is the number of free electrons per cubic meter 
available for conduction, and f is the operating frequency in 
MHz. 

For the vertical case n, = 0, and (4) can be 
utilized to obtain the highest frequency that is reflected by 


the ionosphere at a location when N, is maximun, 


f, = 9(N,, )'” = 


This frequency is called the critical frequency or 
piasma frequency. The critical frequencies of the E, Fl and F2 
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Figure 24. The true reflection height h, and the virtual 
height h'. From [Ref. 11: p. 63]. 


layers are written respectively as f,E, f,Fl and f,F2. 

Many important characteristics of the ionosphere 
are determined from the critical frequencies and the virtual 
heights of the layers. An electronic device called ionosonde 
measures the virtual height h' of the ionosphere by using a 
pulsed-radar technique. This is the height from which the wave 
would appear to be reflected if the ionosphere were a perfect 
reflector. Figure 24 illustrates a ray path through the 
ionosphere and the virtual height h'. If a number of 
measurements are made by the ionosonde, the reflection height 
h against frequency f can be plotted. This is called an 


ionogram. In Figure 25 an ionogram can be seen which indicates 
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Figure 25. A typical ionogram. The superscripts o and x 
refer to ordinary and extraordinary waves. From [Ref.10: 
p- 79}. 

reflection from E, Fl and F2 layers for both ordinary and 
extraordinary waves. 

If the angle of incidence, i,, of the 


electromagnetic wave is not zero, the refractive index, n,, at 


the height at which reflection occurs, is given by 
n, = sin(ip) = [1 - (f, /MUF)? }'? , (6) 


where MUF is called the maximum usable frequency and is the 
highest frequency to be reflected from the ionosphere. 
Frequencies higher than MUF are not reflected, but, instead, 
penetrate the various layers and propagate into space. The MUF 
depends on the specific path, the angle of incidence of the 


shortwave with the ionosphere, the time of day, the season, 
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and the solar activity. The MUF, at any distance, is 


determined by multiplying the critical frequency at zero 
distance by a MUF factor which related to the secant of the 
angle of incidence. Figure 26 shows the MUF factor vs range 
for the various ionospheric layers. The distance that can be 
spanned by the MUF is called the skip distance. 

To account for the MUF variations, the frequency of 
optimum transmission (FOT) is usually used. The FOT is chosen 
to be 85% of the MUF. 

Another quantity is the lowest usable frequency, 
(LUF), which depends on the engineering characteristics of the 
shortwave link. LUF is the frequency limit below which the 


signal-to-noise ratio does not have an acceptable value. 


Ge Variation of the Critical Frequencies 
(i) ___ Solar activity dependence. Sunspots have a 


major effect on the ionization level of the ionospheric 
layers. If the sunspot number is high, the radiation of the 
sun is also high as is the intensity of ionization of D, E, 
and F layers. This means that the attenuation of HF waves in 
the D region increases with the sunspot number, and that 
higher frequencies can be reflected by the E and F regions. 
To measure the sunspot number, the relative 


Zurich sunspot number (Z8N), R, is used. It is given by 


R=k(10g +s) , (7) 


where 
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Figure 26. The MUF factor MUFF, versus range d. The layers 
F2a and F2b represent winter day and night respectively, 
and P2c a summer day or night. From [Ref. 10: p. 82]. 

g is the number of sunspot groups, 

s is the number of single spots, and 

k is the equipment's correction factor. 

The sunspot number fluctuates sharply from day 
to day; thus, for the planning of shortwave links, the twelve- 
month average value Rios is used. Figure 6 shows the dependence 

‘ of the critical frequencies on the average sunspot number R,,. 

Measurements since 1749 show that sunspot 
maxima and minima follow an li-year cycle on the average. 
There is a pronounced correlation of the F2-layer critical 
frequency with the eleven year solar cycle. This can be seen 
in Figure 7 for the monthly mean noon critical frequencies as 


observed at Slough [Ref. 11), over the last solar cycle 
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Figure 27. Variation of the critical frequency fo, with 
the sunspot number. Fron {Ref. 22: p. 23}. 


(1970-1981) . 


{2)___Time variations, The critical frequencies 


vary annually, seasonally and diurnally. Critical frequencies 


the lower Part of the Hr band. On the other hand, the critical 
frequencies of the years of high sunspot activity have much 
higher values. 

The seasonal variations of f£,E and f£.Fl are in 
Phase with solar zenith angle xy. This means that they are 
greater in the summer months. The Opposite takes Place for the 


f.F2 which is in antiphase {13}. This is known as the winter 
@hbomaly and occurs in daytime only. 








| 
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Figure 28. Monthly mean of the noon critical frequencies 
observed at Slough during the last sunspot cycle 1970- 
1981. From (Ref. 11: p. 68]. 


Figure 29 shows diurnal variation of monthly 
median f,F2 for four months and three years of high, medium, 
and low sunspot activity as measured at Slough [11]. Diurnal 
variations are usually observed during winter months. 

9. _ Propagation modes of the Skywave 
The radiation patterns of the transmitting antenna 
and the operating frequency of the shortwave link determine 
which one of the several possible paths the HF ray will 
follow. Figure 30 shows six different ray paths as a function 
of elevation angle for a fixed frequency. The range decreases 
as the elevation angle increases (case 1 to 3) until the skip 


distance is reached. If, however, the elevation angle is 
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Figure 29. Diurnal variation of monthly median f,F2 
observed at Slough. From [Ref. 11: p. 73]. 


further increased, the rays would be trapped in the 
ionosphere, and the attainable ranges would increase rapidly 
(case 4, 5). For still higher angles, a critical angle of 
incidence is reached beyond which penetration of the layer 
occurs (case 6). 

When the shortwave link has several hundred Km 


length, several hops are necessary. The electron 








Figure 30. Ray paths against launch angles for a fixed 
frequency. From [Ref. 11: p. 62]. 


concentration, and thus the ionization of the ionospheric 
layers, is different for the individual sections of the path. 
This means that different transmission conditions take place 
along the path. In Figure 31, examples of possible radiation 
paths for multi-hope links can be seen. 

In multi-path propagation, several signals arrive 
at the receiver with different time spacings. This time 
Gispersion is an inherent characteristic of multi-path 
propagation. The causes of time dispersion include the 
following: 

(a) Ground wave/skywave interference (short links). 
. (b) Skywave reflections from E and F layers. 

(c) Muiti-hop links (see Figure 9). 

(ad) Different elevation angle skywave paths. 


(e) Splitting of the skywave into the ordinary and 


extraordinary wave. 








Figure 31. Various propagation modes for multi-hop links. 
Prom [Ref. 1.2: p. 41). 


Time dispersion is associated with a degraded 
channel performance in terms of error rate. High error rates 
occur as a result of intersymbol interference. 

3. Losses and Field Strengths in a Shortwave Systen 
The received signal power, S in dBW at the output of 


the receiving antenna of a shortwave link is given by 
Se = Pp +m + DD, - Tyee +n, +D, , (8) 
where 
P, = transmitter power (dBW), 


Lot = Path loss in the given direction (dB), 


n, = transmitting antenna efficiency (dB), 
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n, = receiving antenna efficiency (dB), 


D, = transmitting antenna directivity (dB), and 
D, = receiving antenna directivity (dB). 
Assuming that a shortwave system begins with the 


transmitter output and terminates with the receiver input 


(i.e. a closed transmission system), its total loss Lovet in dB 
follows the formula [12], 
Lest = Te, +L, + Uy + ¥p- (G,-G) , (9) 


where 
L,, = losses due to attenuation in free space between 
the transmitting and receiving entennas (dB), 
L, = losses in the ionosphere (dB), 


L, = losses due to the reflection on the ground (dB), 


¥, = fading reserves (dB), 
G, = gain of the transmitting antenna (dBi), and 
G. = gain of the receiving antenna (dBi). 


Referred to isotropic transmitting and receiving 


antennas, the free space attenuation is given by 
Lk, = 20 log(4rd, c/f) , (10) 
where 
d,. = path length in Kn, 


f = operating frequency in MHz, and 


c¢ = speed of light (3x10° m/sec). 








As mentioned above, the emitted electromagnetic wave 
on its way to E or F reflecting layers, passes through the D- 
layer where it suffers substantial losses due to absorption. 
These losses depend on the sunspot number R, and the senith 
angle of the sun x. To calculate the ionospheric losses, the 


absorption index I is used [12], 


I = (1 + 0.0037R) .{cos(0.881x)j'? . (11) 


The absorption index is higher in winter than in 
summer due to the winter anomaly. With the aid of a nomogram 
found in [12] and entering with the absorption index I and the 
elevation angle of radiation, the ionospheric losses L, can be 
determined. 

The reflection of a skywave by the ionosphere changes 
its polarization considerably. Thus, the reflected wave 
reaches the ground with different polarization than the 
emitted one. Assuming that the two components of the energy of 
the incident wave are equal, then the losses due to reflection 


on the ground are given [12] by 


Ly = 10log{(R2 +R) /2) , (12) 


where R? and R,” are the reflection coefficients of the 
vertical and horizontal component respectively. 
The fading of the skywaves is related to changes in 


the ionosphere. Various kinds of fading are defined according 





to their origin: for example, absorption fading, polarization 
fading, etc. Due to fading, the field strength at the 
receiving point is continuously fluctuating. It is found [11], 
that the amplitudes of field strength values undergoing short- 
period changes follow a Rayleigh distribution. To take into 
account the field strength variations, a fading reserve Y, of 
14 GB is generally recommended by CCIR for HF links (see eq. 
9). 

Continuing the analysis of a shortwave system, th. 


received noise in dBW in the receiver bandwidth B is 
N=kT+B+F,+n,+D, , (13) 


where 
kT = thermal noise power density (~-204 dB/Hz), 
B = receiver bandwidth (dBHz), 


F effective antenna noise power factor (dB), and 


D, = receiving antenna directivity factor for the noise 
sources (dB). 
Thus, combining equations (8) and (13), the received 


signal-to-noise ratio in the bandwidth B can be found as 
(S/N)m =S-N=E,-Ly, +D,-N, , (14) 


where 

E, = P,tn,tD, , is the effective radiated power, 
Dr = D.-D, , is the receiving antenna directivity 
factor against noise, and 


N, = kT+B+F, is the noise power at the antenna input. 
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In accordance with the expected grade of service, it 
is necessary to adopt a minimum required signal-to-noise 
ratio. Signals that manage to reach or exceed this threshold 


value are considered successful transmissions. 


C. POLAR AND AURORAL REGION EFFECTS ON HF PROPAGATION 

The high latitude ionosphere is very unpredictable, 
irregular, and variable in nature. These variations and 
irregularities have an important influence on HF waves that 
traverse the auroral and polar regions. In some occasions, the 
HF signals may be attenuated as much as 200 GB, and the effect 
in communications is catastrophic. 

By way of introduction the terms that frequently occur in 
the technical literature will be briefly explained. 

fa) Solar flare: A burst of light occurring in the 
chromosphere of the sun and in proximity to a sunspot. Solar 
flares are accompanied by strong X-ray emission. Flares 
frequently occur at the peak of the sunspot cycle. 

(b) Visual] aurora: The luminosity that can be seen in high 
latitudes (the optical result of the interaction of incoming 
particles with the atmospheric constituents). 

4{¢) Auroral ova]: The region of latitudes containing the 


visual aurora at a specific time. 


(a) Polar cap: The region inside the auroral oval. 





fe) Auroral sgones: The regions of the earth where the 


visible aurora occurs most often; in other words, the auroral 
zone represents the average location of the aurora. 

(f) Solar wind: A stream of plasma that originates from 
the sun. It is composed basically of electrons and protons 
(95%), and alpha-particles (5%). The solar wind is the 
principal medium through which the solar activity is 
transferred to earth; and, as such, is extremely important in 
solar-terrestrial relations. 

{qg) Magnetosphere: The pressure of the solar wind which is 
in motion towards the earth, distorting the terrestrial 
magnetic field (geomagnetic field), and confining it to a 
finite region called magnetosphere. In this region, the 


geomagnetic field exerts a pronounced influence. 


fh) sudden jonospheric disturbances (SID): At times, a 


solar flare produces abnormally high ionization density in the 


D-layer. This phenomenon results in a sudden increase in the 
absorption of shortwaves and fading may occur. These shortwave 
fadeouts (SWF) may last from a few minutes to a few hours and 
occur only on the sunlit side of the earth. In general, SIDs 
depend on the solar zenith angle. SID absorption is the one of 
three types, the others being auroral absorption (AA) and 
polar cap absorption (PCA). 

ona : In addition to the intensified 
ionization produced by a solar flare, electrically charged 


particles cause a disturbance of the earth's magnetic field 
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which may last several days. During these disturbances the 
components of the earth's magnetic field fluctuate over much 
wider limits than under undisturbed conditions. 

There are two types of geomagnetic storms: Those that 
occur sporadically and those that exhibit a strong 27-day 
periodicity. The former type occurs during years of high solar 
activity, but the recurrent type can be found at any phase of 
the solar cycle. Geomagnetic disturbances are more intense in 
the vicinity of auroral zones [13]. 

Geomagnetic activity exhibits a seasonal variation. Figure 
11 illustrates the seasonal distribution of very large 
geomagnetic storms for the period from 1932 to 1984. 
Geomagnetic storms occur most frequently during equinoxial 
months and less frequently during solstitial months. 

(<j) Jonospheric storms: Large geomagnetic storms are 
accompanied by a number of ionospheric phenomena such as: 
reduction of the critical frequency of the Fe-layer (f,F2), 
possible blackout of HF communications and enhancement of the 
F and sporadic E layers. These events are also called 
ionospheric storms. 

(k) Substorm: A substorm is a localized night-time 
phenomenon characterized by the "break up” of visual aurora 
around midnight and accompanied by abrupt changes in radio- 


wave absorption. 
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Figure 32. Seasonal distribution of very large geomagnetic 
storms from 1932 to 1984. From [Ref. 16: p. 5-32]. 


1. Auroral Phenomena 

The polar ionosphere is affected mainly by two 
phenomena: (1) The solar wind induces an electric field in the 
ionosphere which causes ionization drift; as a consequence a 
transport of ionization irregularities through the polar cap 
may occur. (2) Charged corpuscles from the sun and from the 
upper magnetosphere penetrate the atmosphere and cause 
enhanced ionization and irregularities at various ionospheric 
layers. 

The principal high-latitude ionospheric irregularities 


are the following: 








S- Auroral Absorption 


This is the form of absorption which occurs most 


often in polar regions and results from enhanced ionization 
caused when energetic electrons are precipitated into the 
upper atmosphere. The reconnection of magnetic field-lines in 
the magnetosphere is believed to provide such electrons [14] 
which are forced poleward by the earth's magnetic field during 
substorms. 

Auroral absorption is a sporadic phenomenon that 
occurs as a series of discrete absorption enhancements lasting 
from a few minutes to a few hours. It is correlated with 
geomagnetic and auroral activity (though without exact 
correspondence). The auroral absorption has a z2nal occurrence 
pattern with a statistical width of 5° to 8° magnetic latitude 
and is centered somewhere between 64° and 68° of magnetic 
latitude. The precipitated electrons have sufficient energy 
(more than 20 keV) to reach the E and D regions [15]. 

The absorption is at maximum level near midnight 
and in the morning, and at minimum level during the afternoon 
and early evening. The morning absorption appears about an 
hour after a substorm and is not correlated with local aurora. 
On the other hand, midnight absorption is less spread out in 
time and space than morning, and occurs in short-lived but 
intense events which are related to the local substorms [16]. 


Figure 12 gives the probability of a given level of absorption 





at 30 MHz being exceeded at a latitude of maximum absorption 
activity. 

A shortwave passing through an ionized layer 
imparts an oscillatory motion to the electrons which vibrate 
about a mean position. This oscillation causes an increase in 
the number of collisions between the electrons and the heavier 
particles of the ionosphere. As a consequence, energy is 
transferred into the ionosphere and absorbed from the wave. 
The rate of absorption is directly related to the number of 
collisions per oscillation. 

The effect of collisions is to make the refractive 
index n of the medium complex [17]; n* = (¢-jz)?, where p and 
x are the real and imaginary parts of n. For the lower 
atmosphere uw = 1. The total absorption A (in dB), over a path, 
is given [17)} by 


Nv dl 


A= 456x107 :f 
v? + (o + 2, cos6)? 


(15) 
where 

N = electron density in e/m, 

® = angular wave frequency, 

v = electron collision frequency, 


n, = gyrofrequency, and 


® 
tt 


angle of magnetic field with the direction of 


propagation. 
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Figure 33. The probability that a given level of 
absorption will be exceeded at an area of maximum 
absorption. From [Ref. 17: p. 1357]. 


The negative and positive sign in the denominator 
refer to extraordinary and ordinary waves respectively. The 
absorption is thus greater for the extraordinary wave. 

The specific absorption has a maximum at some 
height which depends on the operating frequency. If the 
ionization occurs at great height, v<<w and A is directly 
proportional to (1/0*). If the height is low, v>>w and the 
absorption A is independent of frequency. In practice, 
however, the absorption-frequency dependence is expressed as 


a power law : A(f) = Cf", where n is determined from 


observations with riometers, and C is constant for a given set 
of data. 
b. Polar Cap Absorption (PCA) 

This type of absorption is due to the abnormal 
ionization produced from the arrival of energetic charged 
particles originating at the sun. These charged particles are 
high-energy solar protons (MeV) emitted during solar flares 
and are guided poleward by the earth's magnetic field. Polar 
cap absorption covers the whole region of the polar cap; and, 
unlike auroral absorption, it is smooth and quite uniform. The 
duration of PCA varies from 1 to 10 days but usually is 3 
days. One of the major characteristics of PCA events is that 
they are almost always preceded by a major solar flare. 

The ionization responsible for PCA lies primarily 
in the D-layer of the ionosphere. The effect of PCA on HF 
signals is the very strong attenuation (tens or hundreds dB), 
resulting in the total black-out of communications. 
Statistically, however, the HF circuit reliability does not 
decrease significantly over the ll-year solar cycle. For 
example, during the 1949-1958 solar cycle, the overall outage 
due to PCA was 4 to 5% [16]. Figure 34 illustrates the 
distribution of PCA events for the period 1949 to 1961. The 
association of PCA with the level of solar activity can be 


easily seen. 
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Figure 34. Association of PCA with the level of solar 
activity for the period 1949-1961. From [Ref. 13: p. 279]. 
A remarkable feature of PCA is the marked diurnal 
variation in the magnitude of absorption. Thus, reduction in 
the absorption as much as 4:1 [16] during night hours may 
eccur. This phenomenon is related to the time variation of 
solar proton intensity due to the conditions in the 
interplanetary space [18]. The reduction of absorption 
sometimes makes HF communications possible at night. 
PCA events are also called solar protons events 


(SPE), to emphasize that the solar protons are responsible for 


the absorption. 
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The precipitated soft and hard energetic particles 
in the polar cap and auroral oval regions, are responsible for 
the F-elayer irregularities. There are large and small scale 
irregularities. The large form of irregularity, called a 
"blob", has a horizontal shape and consists of discrete 
patches of enhanced ionization. The electron density of a blob 
may be up to ten times greater than that of the background 
[16]. The size of a blob can reach tens of kilometers. Blobs 
can bend the HF waves from the great-circle bearing by a 
considerable amount. To illustrate this, a blob having a 
density increase of 2x10° cm? in a horizontal range of 300 nm, 
could refract a 10 MHz wave by 40° in bearing at each side of 
the path (80° total) [16]. This represents an increase of 30% 
in path length with respect to the great-circle path; and, for 
a 2000 Km path, it means a multipath spread of 2.0 ms. 

Scatter observations (16] indicate that the blobs 
are usually extended in the East-West direction rather than 
the North-South one. In other words, blobs occasionally show 
variations in latitudinal cross section as a function of 
longitude. 

Instabilities occur at the steep density gradients 
(walls) at the trailing edges of the blob and can scatter the 
HF waves. These are the small scale irregularities with 


subkilometer sizes. 





° = ect 
Sporadic-E (E,) ionization refers to transient 
areas of abnormally intense electron density. These areas are 
horizontally extended and located in the E region. E£, 
ionization occurs primarily during the winter-night period. It 
is very important, because it frequently supports high- 
latitude HF propagation in the 20 MHz region. The auroral E, 
is a significant HF propagation mode for paths up to 2300 Km 
(19}. 
ect sit coughs 
It is known [20] that the inner magnetosphere co- 
rotates with the earth, whereas the outer magnetosphere 
follows circulation patterns determined by the solar wind. The 
difference in the circulation patterns between the inner and 
outer magnetosphere results in the formation of a trough [20]. 
That is, the ionosphere tends to be depleted between the mid- 
latitude and high-latitude regions to form a "belt" of low 
electron density. 
Troughs occur in every season but are a night-time 


Phenomenon at geomagnetic latitudes between 60° and 65°. 
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APPENDIX B ANTENNA RADIATION PATTERNS FOR PROPHET 


‘ AD>777 


ANTENNA LABEL = 777 SLOPING VEE, ANGLE 90 DEG 

TRANSMIT CAPABILITY = Y - POLARIZATION = V BEAMWIDTH = 30.0 
AZIMUTH PATTERN IS UNI-DIRECTIONAL 
deg: -1ll -9 -8 -9-11-13-21-14-13-10 -6-10 -9-26-18-16 -?7 -7-18 -3 
deg: <-8 -6 -7 -7 -6-10-21-25-15 -8 =-4 ~6 6 -8-17-15-14-10 -7 ~4 
deg: <-8 -3 -3 -4 -4 -4 =4 -6-12-14 -9 ~8 3 -1 -1 -3 -7-10 -1 -1 
deg: <-8 -3 <1 -1 -2 -3 ~3 <2 -3 -5 -9-12 “6 -5 -4 -1 1 il-il 1 
deg: -13 -7 -3 -2 0 00221 13 3 3 3 1 ah -2 -4 -5-10 <1 

deg: -22-16-13-10 -8 -7 -7 -6 -6 -5 -4 -3 -2 -l1 -1 0 0 0 0 3 
GW: -38-32-28-26-24-23-22-21-21-20-19-18-17-16-15~15-15-15-15-11~-10 
Freq: 2 4 6 8 10 12 14 16 18 20 30 


AD>778 


ANTENNA LABEL = 776 WHIP ANTENNA, TRANSMITTER 

TRANSMIT CAPABILITY = Y POLARIZATION = V BEAMWIDTH = 360.0 
AZIMUTH PATTERN IS OMNI-DIRECTIONAL 
deg: -40-40-40-40-40-40-40-40-40-40-40~40-40-40-40-40-~40-40-40-40-40 
deg: -40-38-34-31-27-24-22-20-17-17-18-20-21-20-26-19-18-17-18-25-16 
deg: -40-33-29-25-21-18-15-12-10 -8-10~11-13-14-15-14~13-12-10-11-16 
dag: ~40-31-27-22-19-16-13-10 -7 -5 <6 <8 -9-11-l2-12~13-31 -9 -7-21 
deg: ~40-31-27-22-19-15-12 +8 -5 -3 -3 -4 -5 -6 -7 -7 -8 -7 -7 -5 -8 
deg: ~40-36-32-23-24-20-17-13-10 -8 -8 -8 -9-10-10-11~11-10-:0-10 -9 
GWt ~40-40-40-40-37-34=30927-23-21-21-21-22-23-24-24-24-24-23-23-23 
Freq: 2 4 6 8 40 12 14 16 18 20 30 


AD>779 


ANTENNA LABEL = 779 BUTTERNUT (XMITTER) 
TRANSMIT CAPABILITY = Y POLARIZATION = V BEAMWIDTH = 360.0 


e AZIMUTH PATTERN IS OMNI-DIRECTIONAL 
deg: -40~40-40-49-40-40-40-40-40-39-39-39-39-39-40-40-40-40-40-40-35 
deg: -27~24-21-19-18-17-16-15-14-13-12-12-13-13-13-14-14-16-18-15 -9 
& deg: -21~18-16-14-12-11-10 -9 -8 -7 -6 -6 -6§ ~6 -6 -6 -6 -7 ~8-10 -5 


deg: -20~-17-14-13-11-10 -9 -8 -7 -6 -5 -4 -4 ~4 <q -4 -4 -4 <3 -6 <5 
deg: -19~16-14-12-10 -9 -8 -7 -6 -5 -4 -3 +3 ~2 -2 -1 -1 -1 -i 0 0 
deg: -22~20-17-16-15-14-13-12~11-10 -8 -8 -8 ~7 -7 -6 -6 -5 ~-3 -4 -3 

GW: -33~+31-30-28-27-26-25-24~24-23-21-21-21-20-20-19-19-19-18-17-16 
Freq: 2 4 6 8 10 12 14 16 18 20 30 





Figure 35. Radiation patterns, in terms of total gain 
(4B), of sloping-vee (code 777), whip (code 778) and 
"Butternut" (code 779) antennas utilized by Prophet. 
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APPENDIX C NEC DATA SET FOR THE SLOPING-VEE ANTENNA 


CE SLOPING VEE ANTENNA CANGLE VEE 90 DEG.) 
GN1,3,69.62,69.62,-1.,69.62,69.62,0.,.001, 
GH2,3,69.62,69.62,0.,69.62,69.62,.6,.001, 
GH3,41,69.62,69.62,.6,0.,0.,18.,.001 
GNG,41,0.,0.,18.,-69.62,69.62, .6> . 
GW5,3,-69.62,69.62, .6,-69.62,69.62,0.,.001, 
Gens 02 Cena rn Oe rey ek soe re Teo ee tere Se 





Figure 36. NEC dataset for the simulation of sloping-vee 
antenna. 
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